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1 Introduction 



As is known, fermion number may be non-conserved in the standard model at high 
energies, in non-perturbative processes mediated by sphalerons Such processes cor- 
respond to overbarrier transitions between topologically distinct vacuum sectors of the 
theory, where the barrier height is given by the sphaleron mass. As the Chern-Simons 
(CS) number of the gauge field varies, this causes the anomalous nonconservation of 
fermion number. 

In the first quantization picture, the change of fermion number is a result of level 
crossing 0. Consider an adiabatic process, in which the boson fields interpolate between 
distinct vacuum sectors leading to a change of the CS number Ncs- The energy spectrum 
of massless fermions interacting with the boson background then depends parametrically 
on Ncs- When Nqs varies, some of the energy levels cross zero ending up with an 
opposite sign of their energy. If the initial state corresponds to the filled Dirac sea, then 
(anti) fermions will be present in the final state. For the processes involving electroweak 
sphalerons, this picture was studied in Refs.[0],0]. In addition, fermion production in 
dynamical decay of the sphaleron has been considered 

The aim of the present paper is to investigate the level crossing phenomenon arising 
within the context of the Einstein- Yang-Mills (EYM) theory. As is known, in this theory 
there exist regular particle-like solutions 0. Such solutions relate to the critical points 
on the top of the potential barrier dividing distinct vacua in the theory, which allows 
them also to be interpreted as sphalerons One may expect that transition processes 
mediated by these EYM sphalerons become significant at high energies 0. 

It has been shown recently that the EYM sphalerons possess zero energy fermion 



bound state [0. The existence of this state indicates that some of the fermion energy 
levels cross zero when the evolving background field passes through the sphaleron field 
configuration. In the present paper we investigate the whole level crossing picture for 
the case when the background EYM fields interpolate between distinct vacua. 

We define below paths which connect topologically distinct vacua in the EYM con- 
figuration space passing through the sphaleron equilibrium solutions. The field evolving 
adiabatically along such a path is used as the background field in the Dirac equation. To 
discretize the fermion energy spectrum, we compactify the 3-space by lifting the system 
onto a 3-sphere with large but finite radius i?*, and thus the distance between energy lev- 
els is of the order of l/i?*. For this geometrical regularization, we determine analytically 
the free fermion energy levels and zero fermion modes as well. In addition, we derive a 
simple analytic approximation for the whole level crossing spectrum, which gives quali- 
tatively the right result. Finally, we solve the Dirac equation numerically and find the 
spectrum precisely. Our result is depicted in Fig. 3 and shows explicitly the phenomenon 
of level crossing appearing in the EYM theory. 

2 EYM sphalerons and vacuum to vacuum paths 

Consider the action of the EYM theory for the SU{2) gauge group 



EYM 



^ j R^gd'x - i I trF,,F^^^gd^x. (1) 
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Here F^,^ = d^^A^, — d^A^ — ie[A^,A^] is the matrix valued gauge field tensor, e is the 
gauge coupling constant, A^ = A^t'^/2, and (a = 1,2,3) are the Pauh matrices. Our 
sign conventions are those of Landau & Lifshitz . 

We choose the spherically symmetric gauge field as follows 

A = A^dx^ = ^(1 - w)UdU-\ U = exp(z|nV'^), (2) 

where = {sini)cosip, sindsimp, cosd), and the gravitational field is parameterized by 

ds^ = ll { (1 - —)a^de - r\d^^ + smHd^^)] . (3) 

\ r 1 — 2m /r J 



Here le = V'i'Klpi/e is the only dimensional quantity in the problem {Ipi being Planck's 
length); the functions m and a depend on t and r. 

For the static case, the corresponding Lagrangian equations read 

, ,^ 2m, ,2 (w^ - 1)^ , w'"^ 

m' = (1 - —)w'' + \ a' = 2— a. (5) 

These equations are known to possess regular, asymptotically flat solutions discovered 
by Bartnik and McKinnon 0. These BK solutions are labeled by an integer, n. The 
function w{r) for the ra-th solution, Wn{r), has n nodes in the domain < r < oo and 
satisfies the following boundary conditions: 

w„(r) = 1 + O(r^) as r ^ 0, w„(r) = (-1)" + a„/r + 0(l/r^) as r ^ oo, (6) 

where a„ are numerically known constants |0. The corresponding metric functions 
m„(r) and cr„(r) increase monotonically from m„(0) = to m„(oo) = and from 
<7n(0) = cr„ to (Tnioo) = 1, respectively (see Fig.l). The ADM masses of these solutions 
are M„ = A/47rMp;m„/e, with Mpi being Planck's mass, where m„ increases as n grows 
from the minimal value mi = 0.828 to m^ = 1. 

As was shown in Ref. the odd-n BK solutions may be treated as sphalerons. 
This interpretation is based on the fact that, through any such solution, one may find 
a one-parameter family of fields which interpolates between topologically distinct EYM 
vacua. If A denotes the parameter of the family, the corresponding gauge field can be 
chosen as follows 

A^dx^ = —(1 - Wn{r))UdU-\ U = exp(«AnV'^), (7) 

where A ranges from zero to it. The metric functions m and a in (|^) are given by 

a{r) = exp{—2sin^ X j w^ — }, 

Jr T 

m{r) = ^ + szn'X^^^^i—^)adr. (8) 

a[r) Jo 2r^ 

It is worth noting that, for any value of A, the fields of this family satisfy the two Einstein 
equations Gq = SttGTq and G^. = SttGT^; when X = 7i/2 these equations reduce to those 



in (|^) 0. Obviously, for A = 7r/2, the fields d^);® coincide with the n-th BK solution 
field. When A approaches zero or vr values, the spacetime metric becomes flat and the 
gauge field (|^ vanishes. Thus, the family of fields (^,(||) describes a loop (in the EYM 
configuration function space) which interpolates between the trivial vacuum and the n-th 
BK solution when A runs from to 7r/2, and returns back to the vacuum as A changes 
from 7r/2 to tt. 

It is important that this loop can be transformed to a path interpolating between 
neighbouring EYM vacua with different winding numbers of the gauge field. The gauge 
transformation 



A^^UiA^ + -d^)U-\ with [/ = exp(z^K - IKO (9) 



yields 



A^dx^ = -(1 - Wn{r))U+ dU-' + -(1 + Wn{r))U^ dUz\ 

U± = exp{iX{wn ± l)nV72), (10) 

whereas the gauge invariant metric functions (§) remain unchanged. One can see that, 
when A runs from zero to tt, the field (0) interpolates between two different pure gauges. 
We will refer to the gauge field potentials (|^) and (|10D as gauge field represented in the 
gauge / and in the gauge // , respectively. Notice that, in the gauge // , one has 
Aa = 0{r) when r — > and Aa = 0(l/r^) as r — >■ oo, with Aa being the field component 
with respect to an orthonormal frame. 

Let A depend adiabatically on time, in such a way that A(t = —oo) = and 
A(t = oo) = 71. Then, in the gauge /, the time component of the field appears, whereas 
in the gauge // one has Aq = 0. Recall the Chern-Simons (CS) current and its divergence 



2 ^fiua/s 2ie ~ 

K^ = —tr^A,{VaA(,-—A^Ap), V^K^ = -—rtrF^^F^' , (11) 



where Vq is the covariant derivative, and F^^ is dual tensor. Integration leads to the 
following relation: 



e2 rt 



167r2 



dt I d^x^/^trF.^Ff"' = I d^Xy/^K° * /"* - /" 



+ / dtfKdJ:. (12) 

oo J 



In the temporal gauge II , the surface integral entering the right hand side vanishes, 
while the gauge invariant left hand side is M 



'6 /■* • f°° 3 1 

— / dt Xsin^X / dr w'^{w'^ — 1) = — (A — sinXcosX){-wf^ — w. 
2tx J-oo Jo Atx 3 



ui„(oo)=-l 

■ (13) 

Wn(0) = l 



Noting that the gauge potential (|TOD vanishes at A = 0, one obtains the following value 
of the CS number along the path (|10D §, 1: 

- —(A — sinXcosX), (14) 

TT 



Ncs = J V^K'd'x 



which changes from zero to one as A runs from to vr (see also Ref.[|TI|] for a quite 
different approach to the calculation of the CS numbers). 
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Thus the fields (||), (|T0|) interpolate between two EYM vacua with different 
winding numbers of the gauge field, and coincide with the (gauge transformed) field of 
the n-th BK solution when A = 7r/2. This allows one to treat odd-n BK solutions as 
sphalerons, "lying" on the top of the potential barrier separating distinct topological 
vacuum sectors of the theory. The profile of this barrier may be obtained from Eq.@ 
because the metric function m{r) obeys the initial value constraint, which allows us to 
define the ADM energy |Q: 

f/(A) = lim m(r) = 



■>oo 

sm A 



/■oo 7/) 1"^ (IT 

. / iw'^ + WA ^ " , ' )exp(-2stn^X / w'^—)dr, (15) 
JO 2r^ Jr r 

(a plot of the function U{NcsW) for = 1 is depicted in Fig.2). 



3 Compactification of the 3-space 

Our aim is to investigate the energy spectrum of chiral fermions in the external fields 
(0)' (0)' ® which form the potential barrier (^). As one expects, the change of the CS 
number across the barrier causes a non-trivial spectral flow. This concept is, however, 
not meaningful since the energy spectrum is purely continuous. Therefore, we apply 
some regularization procedure to discretize the spectrum. The most natural way is to 
compactify the 3-space by lifting the system onto a 3-sphere with large but finite radius 
i?*; the limit R^: ^ oo being eventually implied (see also Ref. for other choice of the 
regularization scheme). 

For this we multiply the spatial part of the line element (||) by an appropriate 
conformal factor 

ds^ = ll [{I - '^)aV - (l + ^) " { + r'^d^' + --^^V)}) , (16) 

with functions m and cr still given by Eqs.(§). When parameter A in (^ interpolates 
between zero and vr values, this line element describes an everywhere regular geometry 
on a manifold whose spatial section is topologically S^. When A = 0, tt, one has m = 
0, cr = 1, so that the spatial section is the 3-sphere with radius and standard metric. 

The lift of the gauge field is locally well defined by Eqs.(^ or (0). However, in 
the gauge / , the field is singular at f = 0, where f = 1/r (we will call this point north 
pole). Consider the gauge II . In the vicinity of the north pole, the components of the 
field ( p!OD with respect to an orthonormal frame are 

Aa = ^aa{^,^) + Oif), (17) 

with tta being some functions of the angular variables, the constant is the same as that 
entering Eq.(|^). One can see that this field is not defined uniquely at f = 0. One can not 
improve the situation just by performing a further gauge transformation in the vicinity 
of the pole (because the maximal decay rate of the field (^^), Aa ~ 1/?"^, is fixed by the 
asymptotical behaviour of the corresponding gauge invariant energy density: Tq ~ 1/r^). 
Notice, however, that for our geometrical regularization, the normalized fermion modes 



4 



always vanish at the north pole (see Eg. (^7]) ) , so that this ambiguity does not enter the 
Dirac equation. 

One can understand this also as follows. Let us introduce for a moment an addi- 
tional regularization, demanding for the amplitude Wn, in the vicinity of the north pole, 
to tend faster to its asymptotical value: 

Wn{r) = (-1)" + a„/r^+^ + ©(l/r^), where e ^ +0. (18) 

The gauge field is then everywhere regular, whereas the field topology, that is, the CS 
number, remains the same. It turns out that, in the vicinity of the pole, the regular 
power series solution for fermions converges uniformly as e to that corresponding to 
e = 0. Thus fermions do not feel the e-dependence of the field. 

We assume therefore that Eq.(p!0|) specifies the regular gauge field on a compact 
manifold whose geometry is described by the line element (|16D. 



4 Fermions on an EYM background 

Our next task is to perform separation of variables in the Dirac equation. For the 
sake of generality, we will consider in this section massive doublet fermions interacting 
with arbitrary time dependent, spherically symmetric EYM background fields. Some 
particular cases of this problem were considered previously in Refs. [l^ , |]10 . 



Consider the action of an SU{2) doublet of fermions in an EYM background 
Sf = + eA,,) - M) ^y^d^x, 



(19) 



with M being fermion mass; this action is invariant under the gauge transformations 



A^-^U{A, + -d^)U~\ ^ 



U^, U E SU{2). 



(20) 



It is most economic for our purposes to make use of the Newman-Penrose (NP) 



formalism of spin coefficients |T2[. This formalism is known to be a convenient tool for 
studying the wave equations of fields of various spin in the algebraically distinguished 
spaces of general relativity [Q, |TB|. In the spherically symmetric case, the role of the 
spinor angular functions is played then by the spin- weighted spherical harmonics. (Notice 
that, in our problem, the separation of variables in terms of spherical spinors can also 
be performed |]T6[, but the calculations are somewhat tedious.) 

In spinor representation, bispinor ^ decomposes as = r .In the chiral 

limit, two component spinors and ^a' (they carry also isospin indices) become right- 
handed and left-handed parts of bispinor respectively. The Dirac equation reads 



(V 



AA' 



ieAAA')<P' 



leA^^'YU' = 0^ (21) 



where the spin-tensors Vaa' and Aaa' are related to the covariant derivative operator 
and the gauge field potential A^, via van der Waerden connection quantities a^^, 

V„, V/i 4' = a'i .,V„,, A„, ^ A 



AA' 



'-'AA'^fJ.- 
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The conserved current is 

which gives the inner product 

< vl/^lvl/i >= J J^^gd'x = I a%,{cl>U2 + iM')V^9d'x. (22) 

Introduce a spinor dyad {o^^l^), oal^ = 1, which defines a null NP tetrad 
{l,n,m,m) according to the scheme 

P^o^o^', n'^^A^', m'^^oV, m^^6^o^', 

non-zero scalar products being l^rif^ = —m^fh^ = 1. Expanding (j)^ and ^a' over the 
spinor basis, 

0^ = $0 + $1 U' = So OA' + Si LA', 

and using the NP quantities 

o^'Vaa'O^ = e-p, o^'Vaa'^^ = vr - a, i^'V^A'O^ = /? - r, i^' Vaa'^^ = /i - 7, 
one obtains the Dirac equation as 

M 

{D + e- p - leAi) % + (5 + ti - a - ieArn) $i = «^ Si, 

v2 

M 

((5 + /3 - r - ?eA^) $0 + (A + - 7 - ^eA„) $i = -z-^ Sq, 

M 

(D + - p* - «eA,) So + (5 + TT* - a* - leAm) Si = $1, 

v2 

M 

(5 + /3* - r* - zeA^) So + (A + /i* - 7* - leA^) Si = z-^ $0. (23) 

Here the standard definitions used are: D = l^d^, A = n^d^, 6 = m^d^, and also 
Ai = I'^A^, An = n^A^, Am = m^A^, Am = ffi^A^ = A]^] asterisk denotes complex 
conjugation. 

Let us choose the spacetime metric in the form 

ds^ = e^^dt" - e^^dr^ - R^{dd^ + sin^M^^), (24) 

where functions a, h and R depend on t and r. The general spherical gauge field fl^ can 
be represented as follows: 

eA = WqLi dt + WiLi dr + {p2 L2-(l-pi) L3} di3 + {{l-pi) L2+P2 L^}sin^ dip, (25) 

where 

Li = -n"'T°' = sindcosip Ti + sindsiriLp T2 + cos'd T3, L2 = d^Li, L3 = -d^pLi, (26) 

2 sinv 

with Ta = t'^/2; Wq, Wi,pi,p2 being functions of t and r. The gauge transformation ( PU] ) 
generated by 

U = exp{zn{t,r)Li) (27) 
6 



preserves the form of the field p5|), altering the functions Wq, Wi, pi, p2 as 

Wo^Wo + Cl, Wi^Wi + n', p± = pi ± ip2 exp{±in)p±; (28) 

here dot and prime denote differentiation with respect to t and r, respectively. 
We choose the null NP tetrad as 

F = (e- e-^0,0), n'^ = i(e-^-e-^0,0), = -1^(0, 0, 1, -^), (29) 
the non-vanishing NP coefficients being 

p=-^{Re-'' + R'e-'), f, = ^{Re- - R' e-') , ^ = -l(6e-'^ - aV), 

e = l{be"^ + a'e-'), = (30) 

The tetrad components of the gauge field potential are 

eAi = {e-''Wo + e-'W,)L3, eA„ = ^{e-''Wo-e-'Wi)L,, leA^ = -l-(p+-l)L+, (31) 
where L± = Li ± zI/2. Using Eqs.(|5IlD,(PT|), the Dirac equation ( PHj ) takes the form 
^+'^'0 + ^(A- + ^ctg^ + (p- - 1)L_) $1 = Si, 
/2 1 

V-^i + ^(A+ + -ctgi^ - {p+ - 1)L+) $0 = -tV2M Sq, 

P+So + ^(A+ + ^ct(7i9 - (p+ - 1)L+) Si = $i, 
\/2 1 

P-Si + ^(A_ + -ctg^ + (p_ - 1)L_) So = %V2M $0, (32) 



where 



R ' 2 



A± = de ± -^9^. (33) 



Introduce a pair of isotopic spinors: 





/ 
















V 







-sinf e-'^/2 



(34) 



with the following properties 

A±x^ = \ctgdx^, (A± + ^ct(77?)x^ = TX^. (35) 



To decouple angular variables in Eqs.(|32D we make the ansatz 
1 



R 



(36) 

where sYjm{'&^^) are spin weighted spherical harmonics ||l3l, and the functions F^, 
[a = 0, 1) depend on t and r alone. A gauge transformation of the form (pOD,(p7D induces 
the replacement 

F± ^ exp{±zn/2)F^, Q± ^ exp{±in/2)Qf. (37) 
Taking into account Eg s . (1351) and also the following standard relations for the functions 



{d{^ =F —^^'P =t -5 C^fi'^) sYjm = ±\l {j ± S)(j =F S + 1) sTlYjm, 

one can see that the ansatz (^) solves Eqs.(|32D, provided that the following equations 
hold 

DlF+ + -F^ = iMQt, 



D_ 






KFo 


+ R^^ 




D*_F+ 














+ ^Wo 





-iMF^, 



' 



R 



; 



(38) 



where A = wj(j + 1), and 



D± = e-'^(9, + ^ + Vo) ± e-^(9, + Vi). 

For the S'-modes, j = 0, the spherical harmonics ±iYjm vanish, and the resulting equa- 
tions reduce to the Fq' = Fi = Qf = Qq = A = case. One can see that the 
transformations (PHD,(P7D leave Eqs.(P^) invariant. 
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5 Fermion energy spectrum 



For the background field quantities entering the system (|38|) , we choose now (|]) , (|ToD and 
([TE|), describing EYM fields of a vacuum-to- vacuum path. We treat these fields, in the 
adiabatic approximation, as a sequence of instantaneously static configurations, with A 
being a parameter. Our aim is to study the energy spectrum of fermions as a function 
of A. 

In the chiral limit, the two chiralities decouple. In the following, we consider just 
the right-handed spinor component. We are interested in the 5- wave sector dynamics, 
because, as we expect, it is in this sector where the level crossing phenomenon will occur. 
The following combinations are introduced for the fermions functions surviving in the 
case: 

+ F+ = iV2e-"^'ij+{r), - F+ = V2e-'"^'ip- (r). (39) 
The non-trivial equations (^) take now the form 

where = — 2mr; the functions m and a are given by (|^). The quantities specifying 
the background gauge field ( |I0| ) (taken in the regular gauge II ) read 

Wi = Xw'^, pi = cosXcosWnX + WnSinXsinWnX, P2 = cosXsinWnX — WnSinXcosWnX. 

(41) 

Due to (|28|) , (|37|) , these equations possess the following gauge symmetry 
ip ^ ip cos— — ip sin—, If) tp sin— + ip cos—, 

Wi Wi + fl', pi pi cosQ — p2 sinQ, p2 pi sinQ + p2 cosQ, (42) 

where ^2 = Q{r). Sometimes it is convenient to utilize this symmetry in order to search 
for the solutions to Eqs.(^OD, working not in the regular gauge // for the background 
gauge field, but in the gauge / instead. The corresponding field parameters are simpler 
in this case: 

Wi =0, Pi = cos'^X + WnSin'^X, p2 = {I — Wn)sinXcosX . (43) 

Once the solution is found, one has to switch back to the regular gauge // using the 
transformation (|4^) , with the parameter 

n = X{wn - 1). (44) 

The conserved inner product (0) reads 

2 / 2 \ ~^ 

< ^2|V'i >= /^^ ^ (1 + ^) {^t^t + '^2i^i)- (45) 
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The formal series expansions around r = and r = oo give the following boundary 
conditions for the normalizable solutions to Eqs.(^), (^) 

iP+ = 0{r^), ip- = 0(r), as r 0, (46) 

1 

^+ = 0{-), ^- = 0(^1, as r^oo. (47) 



For each value of A, the solutions with the asymptotics ( ^61) and (|47|) have to be matched 
by adjusting the value of u. As a result, one obtains the energy spectrum, uj{\). 

Before passing to tackling the problem numerically, it is helpful to investigate some 
particular cases. Firstly, consider the beginning of the vacuum-to- vacuum path, A = 0. 
Then, both in the / and // gauges, one has pi = 1, p2 = Wi = 0, that is, the gauge 
field is zero; also m = 0, a = 1. Eqs.(^OD describe in the case free fermions on i? x S^. 
Introducing a new radial coordinate x G [0, tt] by r = i?^, (1 + r^/4/2^) sin x, one can 
represent Eqs. (^OD as 

1 \ , u ( d l\,z/ 

- -^^ = 0, : ^+ + -i)- = 0, (48) 




sin X J 2 \dx sin x J 2 

where u = 2ujR^. The general solution of these equations reads 



A 



I vcos{vx/2 + B) - ctg{x/2) sin{iyx/2 + B) 
\ -vsin{vx/2 + B)+ tg{x/2) cos{ux/2 + B) 



(49) 



with A and B being integration constants; i/ 7^ ±1. This solution vanishes at x = only 
if i? = 0. When v is an odd integer, the solution vanishes also at a; = tt. The values 
z/ = ±1 correspond to the degenerate case. The solution then includes power terms, and 
it can not be regular both at x = and x = n. One obtains therefore the following 
spectrum of fermions (see also |llO| ) 



^(A = 0) = -^, z/ = ±3, ±5.... (50) 

Consider now the end point of the vacuum-to-vacuum path, A = tt. To find the 
spectrum, we observe that, in the gauge / , the gauge field parameters (|43| ) take the 
same values both at A = and A = tt. This means that the solution ( P9| ) can be used 
also when A = vr, but only in the gauge / . To pass to the gauge // , one just performs 
the transformation (|4^), so that the spectrum (|50|) remains the same. 



We conclude therefore that, as A varies from zero to vr, the spectrum ([50| ) maps 
into itself. Certainly, this does not mean that each energy level maps into itself. The 
level number, z/, may change. For instance, some of the levels may cross zero, ending up 
with an opposite sign of their energy as A ^ tt. Observe that, in the gauge / , Eqs.(^D|) 
are manifestly invariant under 

oj —u, A ^ vr — A, ip^ — > —ip^, ip~ (51) 

This means that normalizable solution with zero energy may appear only if A = 7r/2, 
and one of the functions tp^ vanishes. In the gauge / , the corresponding solution (first 
discovered in Ref. ||10||) can be found from (^): 



= 0, ip = exp [ j ) ; (52) 
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to pass to the regular gauge // , one uses (0), 

Thus, we can see that the spectrum maps into itself as A runs from zero to vr, 
and also some of the levels cross zero at A = 7r/2. To get an idea about the whole 
level crossing picture, we use the following simple analytic approximation which gives 
qualitatively the right result. Let us put in P0|) for any A, m = 0, a = 1, that is, we 
neglect the change of the geometry along the vacuum-to-vacuum path. Also, we replace 
the function Wn in (p!]), ( ^31) by a step function w obeying the same boundary conditions: 



Wnir) w{r) = 1, r < tq; w{r) = —1, r > tq. (53) 
Notice that the field topology remains intact, as the Chern-Simons number still changes 



in accordance with Eq.(n), for this number depends only on the boundary conditions 
for Wn (see Eq. ([T3|) ). The normalizable zero energy solution ( p2|) also exists in this case. 

To find the solution for the fermion wave function, we observe that, in the gauge 
/, the function w describes a field which vanishes when r < tq and is a pure gauge 
when r > To (see Eq.(|^). The resulting solution is therefore obtained by matching the 
free fermion modes in the domain r < tq and the gauge transformed free modes in the 
region r > tq. Let us use the same independent variable x = x{r) as in (^). When 
X < xq = Xo(ro), the regular solution is given by Eq.(^) with B = and A = 1. To 
obtain the solution in the region x > xq, one takes the solution (^91) with arbitrary A 



and B, and performs the gauge rotation (52) with Q = 2X (to find Q one compares (^ 



(|53| ) and (|27|) ). Matching these solutions at x = xq specifies the values of A and B. 
Next, the solution in the region x > Xq has to be rotated to the gauge // and required 
to vanish at x = vr, which gives the quantization condition. It suffices to demand for 
only one component to vanish, say as the other will then vanish automatically, due 
to (^71) . In the limit xq — > 0, the quantization condition takes the following simple form 

ctgX = ^z/xo(l + 0{xo)) + -^xl{u^ - 1)(1 + 0{x^)) tg{u^). (54) 

This equation gives the spectrum a;(A) = v{X)/2R^. When A = 0, vr, ctgX diverges, which 
requires v = ±3, ±5,..., in agreement with Eq.(^). As A runs from zero to vr, ctgX 
decreases from +oo to — cxd, which forces z/(A) to vary from u{0) 7^ 3 to //(vr) = u{0) — 2, 
and from z/(0) = 3 to /^(vr) = —3. 

We arrive, therefore, at the following level crossing picture for the right-handed 
chiral fermions. When the CS number varies from zero to one, each energy level shifts 
downwards replacing the preceding one with the lower energy. The lowest positive energy 
level crosses zero when Ncs = 1/2 and becomes eventually the highest negative energy 
level. If the initial state corresponds to the fermion vacuum, then one (left-handed) 
antifermion appears in the final state. 

Finally, we turn to the numerical analysis. We consider the vacuum-to- vacuum path 
([T0|), ( p!6D passing through the (lifted) ground state EYM sphaleron field configuration. 
The compactification radius is chosen to be = 10^ which is much larger than the 
sphaleron size (the latter is of the order of unity; see Fig.l). Numerical integration of 
the Dirac equation (|40D with the boundary conditions (|^), (^Tf ) gives the spectrum. The 
result is depicted in Fig. 3 and confirms the qualitative picture described above. For the 
paths passing through the higher (n > 1) EYM sphalerons, the level crossing picture is 
qualitatively the same. 
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6 Conclusion 



Our results show exphcitly the anomalous non-conservation of fermion number in the 
transition processes mediated by the EYM sphalerons. It was argued recently that such 
processes can be significant at high energies due to interference effects 0. One may 
expect then an efficient fermion number non-conservation in the high energy limit. To 
determine more precisely the production of fermions, by going beyond the adiabatic 
approximation, requires a numerical solution of the full time dependent problem 

Our results may be of use also for further investigation of the quantum correction 



to the classical energy of the sphaleron barrier, arising from the fermion see 0] , [O] . 
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Figure captions 



Fig.l. Behaviour of the gauge field function, Wi{r), the metric functions, mi(r), o"i(r), 
and the radial energy density r^Tg for the ground state (n = 1) EYM sphaleron solution. 
In the units used, the sphaleron has the typical size of the order of unity. 

Fig. 2. The potential barrier dividing neighbouring topological EYM vacua. The top of 
the barrier is the position of the ground state EYM sphaleron. 

Fig. 3. The fermion energy spectrum along the path which interpolates between distinct 
topological EYM vacua passing through the ground state sphaleron field. 
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